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The transition from a gapless liquid to a gapped dimerized ground state that oc- 
curs in the frustrated antiferromagnetic Majumdar-Ghosh (or J\ — J2 Heisenberg) 
model is revisited from the point of view of entanglement. We study the evolution 
of entanglement spectra, a "projected subspace" block entropy, and concurrence in 
the Schmidt vectors through the transition. The standard tool of Schmidt decompo- 
sition along with the existence of the unique MG point where the ground states are 
degenerate and known exactly, suggests the projection into two orthogonal subspaces 
that is useful even away from this point. Of these, one is a dominant five dimensional 
subspace containing the complete state at the MG point and the other contributes 
marginally, albeit with increasing weight as the number of spins is increased. We find 
that the marginally contributing subspace has a minimum von Neumann entropy in 
the vicinity of the dimerization transition. Entanglement content between pairs of 
spins in the Schmidt vectors, studied via concurrence, shows that those belonging to 
the dominant five dimensional subspace display a clear progress towards dimeriza- 
tion, with the concurrence vanishing on odd/even sublattices, again in the vicinity 
of the dimerization, and maximizing in the even/odd sublattices at the MG point. 
In contrast, study of the Schmidt vectors in the marginally contributing subspace, 
as well as in the projection of the ground state in this space, display pair concur- 
rence which decrease on both the sublattices as the MG point is approached. The 
robustness of these observations indicate their possible usefulness in the study of 
models that have similar transitions, and have hitherto been difficult to study using 
standard entanglement signatures. 
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I. INTRODUCTION 



Entanglement in many-body systems has been extensively studied recently [iMa] ever 



since the remarkable properties of quantum entanglement have come to^ 
especially through its various uses in quantum information processing 



phase transitions 



to 
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3e understood, 



llj . Quantum 



12| which occur at zero temperature as some external parameter is changed 



has been particularly addressed with the help of entanglement 13Hl6|. The Ising model 
critical point for example has been shown to have an entropy of entanglement that scales 
logarithmically (~ In L) with the length, L, of the spin chain while away from criticality it is 



independent of L 17H19l| . While many condensed matter systems have been studied with the 



help of such a "block" entanglement entropy 
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19M23H. it seems more natural to consider 
2m in contexts where dimerization 



measures of two-body entanglement like concurrence 
occurs One of the well-studied Hamiltonians in this context is the Majumdar- 

Ghosh, or the J\ — J 2 Heisenberg model 29, 3o|. This model has a well known transition 



rom a gapless critical phase to a gapped phase with short range correlations and dimer order 



31 



32] . Earlier studies on this model from the entanglement perspective have employed 



scaling behaviors of the von-Neumann entropy of contiguous blocks of spins, the valence- 



33|, 



34j, and other measures of multipartite entanglement 



band entanglement entropy 15, 
351 ] to study this transition. 

Direct signatures of the dimerization transition in this model using two-spin entanglement 
measures such as concurrence have been elusive. 

In the J\ — J2 model the antiferromagnetic nearest neighbour Heisenberg chain is aug- 
mented with a next-nearest neighbour Heisenberg interation which is also antiferromagnetic 
29, ai 3. The interest in this has been considerable since Majumdar and Ghosh pro- 
posed this as a model with an exactly solvable ground state that shows dimerization at 
J2/J1 — 1/2, the so-called MG point. The Hamiltonian is 



N 



N 



H — J\ 2J &i ■ Pi+l + J2 2J <?i ' ^+2> 



(1) 



i=l 



i=l 



with Ji, J2 > 0. The particles are of spin- 1/2 and <t, are the Pauli matrices. Periodic 
boundary conditions are assumed, so that <tn+i = o\ and c?at +2 = 0%. At the MG point the 
ground state is doubly degenerate and ground state manifold is spanned by two states \Rn) 



3 



and \L N ), where 

\R N ) = (1 2) (3 4) • • • (N - 1 N), \L N ) = (2 3)(4 5) ■ ■ • (N 1). (2) 

Here, for example, (12) refers to the singlet state ^=(|01) — |10)) of spins 1 and 2, and 
|0) and |1) are eigenstates of a z with eigenvalues 1 and —1 respectively. Thus at the MG 
point the degenerate ground states can be considered to have maximal nearest neighbour 
entanglements as the entanglement of a singlet is the maximum possible between two spin- 



1/2 particles. Thus it would seem natural that entang 



ement between spins is enhanced at 



a transition from a spin-liquid to a dimerized phase [31|, |32|. a transition which occurs when 
J 2 /Ji ~ 0.24.... While a fair amount of literature already addresses this 15, [35I H, the 
present Paper revisits the issues 

from the point of view of analysis of entanglement spectra at finite lattice sizes, a non- 
standard block entropy, and concurrences in the Schmidt vectors of the reduced density 
matrices. We expect this to be of interest in larger classes of problems where the possibility 
of transitions from spin liquid to dimer order need to be investigated. 

The strategy is to focus on the fact that this happens to be one of the rare systems where 
at least at one point in the phase diagram, namely the MG point, the ground state can be 
solved exactly and has a form simple enough to enable the evaluation of the entanglement 
spectrum analytically. The nature of the entanglement spectrum at this point suggests a 
separation of the state into two orthogonal components with supports in what one may call 
a "MG" and a "non-MG" subspace. At the MG point, the MG subspace solely contributes 
towards the construction of the ground states, hence the terminology. The MG subspace is 
only 5-dimensional and in the range < J 2 < 1 /2 seems dominant at least for small system 
sizes. In fact the competition between this subspace and its complement seems to be crucial 
for the emergence of a dimer order. For later convenience the non-MG subspace is denoted 



as MG. However these subspaces are not unique in a way that is elaborated in the next 
section. 



It is found that a suitably defined entanglement corresponding to the MG component 
of the wave function has a minimum in the vicinity of the dimerization transition. Thus 
while there do not seem to be simple signatures (except for scaling with L) in either the 
entanglement entropy of the state 34j, or its dominant part, namely the MG component, the 



typically small MG components apparently carry information that may signal the transition. 
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The separation of the entanglement spectrum into these two components also allows for a 
detailed study of the entanglement of the eigenvectors of the reduced density matrices. While 
much attention has concentrated on the entanglement spectrum per se, it is but natural that 
the eigenvectors have significant information in them. The entanglement in these Schmidt 
vectors is studied, especially the concurrence between nearest neighbors. It is observed 
that for vectors corresponding to the MG component, a clear dimerization happens, with 
alternate pairs of nearest neighbour entanglements either increasing to the maximum value 
as the MG point is approached, or vanish in the vicinity of the dimerization transition. 



II. THE EIGENVALUE SPECTRUM OF THE GROUND STATE OF THE MG 

MODEL 

First, a separation of the ground state of the MG Hamiltonian, say \^(J 2 )), (Ji = 1 from 
now) into two distinct orthogonal states, with properties described below is sought. Thus, 

|*(J 2 )) = a(J 2 M MG (J 2 )) + /3(J 2 )fe(J 2 )>, (3) 

where IV'mg^)) approaches a superpostion of \R) and \L) as J 2 — > 1/2, and a(J 2 ) — > 1, 
P{J"i) — ► in the same limit. The state IV'mg^)) is orthogonal to this, and will play a 
rather important role here. This non-MG part forms a small fraction of the whole state, 
at least for small N (e.g. 3% for N = 16) . With increasing number of spins though this 
component grows and the detailed manner in which this happens as a function of J 2 is 
interesting and may hold information about the dimerization. However, by definition this 
component decreases to zero at the MG point of J 2 = 1/2 for all N. Such a separation 
is possible, but is potentially non-unique, as demonstrated further below. Throughout this 
paper the number of spins is an even number, and there are two main subclasses: N/2 even 
and N/2 odd which are simply referred to as "even" and "odd" cases. Also from the point of 



view of symmetry, the translation symmetry is broken in the projected MG and MG parts 
for all J 2 . 

One would especially like to treat the interval < J 2 < 1/2, which contains the point 
where there is a gapless to gapped transition. When J 2 = the ground state (and, indeed 



any, excited state) is solvable via the Bethe ansatz 38), however the explicit forms are 



unwieldy and difficult to analyse in detail. Thus a rather "complex" antiferromagnetic 
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ground state at J2 = evolves to a rather simple dimerized state at J2 = 1/2. That a part 
of the ground state can be identified for all Ji in the interval [0, 1/2) that evolves to the 
dimers at the MG point is not necessarily obvious and is elucidated in this paper. 

That this is possible is strongly suggested from a study of the Schmidt decomposition 
of the ground state. Let iV spins in the chain be split into two parts (say A and B) of 
contiguous spins having Na and Nb particles each. This paper will concentrate on the 
cases when Na and Nb are even numbers as well. This ensures that the subsystems under 
consideration are of the same parity (number of spins odd/even) as the original chain. The 
Schmidt decomposition in terms of vectors from these two halves reads: 



2 N A 



= E ^W2)\W2))a\W2))b. (4) 
3=1 

Here Aj(J 2 ) are the eigenvalues of the reduced density matrix (RDM) 

^ A (J 2 )=tr B (|*(J 2 )>(*(J a )|), 
and \4>j{J2))A are the corresponding eigenvectors. The eigenvalues Aj( J2) are also dependent 



on the partition size Na, but this is not explicitly indicated. The von Neumann entropy 
Sn a (J%) = — Sj=i Aj(^2) log(Aj( J2)) is a measure of the entanglement between par 
and B. There have been several works that study the so-called entanglement spectrum 



39- 



4lj which is defined as {— ln(Aj), j = 1,2,...} in many systems, such a spectrum naturally 
containing much more information than just the entropy. For most part of this paper, unless 
otherwise mentioned, Na = N/2 for the even case and Na = N/2 — 1 for the odd case. It 
must be noted that only for the entanglement spectrum we take logarithm to the base e and 
for all other measures of entropy the logarithm is taken with respect to base 2. 

Fig. [1] (top row) shows the eigenvalues of the RDM, where the principal eigenvalues 
corresponding to the MG subspace are seen clearly. The largest eigenvalue decreases as 
the MG point is approached from the Heisenberg. The second largest eigenvalue actually 
comprises of a triplet that is almost a constant as J 2 varies in [0, 1/2). The smallest of the 
eigenvalues that is clearly visible in this figure increases as the MG point is approached and 
indeed seems to become significant in the vicinity of the dimerization transition. It is shown 
below that at the MG point this eigenvalue is coupled with the largest one. As N — > 00, it 
approaches the value 1/8. The Schmidt vectors (pure states of N/2 particles) corresponding 
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FIG. 1. The eigenvalues of the RDM pat a (J 2 ) for N = 16, 24 (corresponding to N A = N/2) and 
N = 18, N = 22 (corresponding to Na = N/2 — 1) are shown. Top row: 50 largest eigenvalues 
are plotted. Prominently seen are the 5 "dimer" or MG subspace eigenvalues, a large eigenvalue 
around 0.6, the almost constant triplet around 0.1, and the small but rising eigenvalue that becomes 
important around the dimerization transition. Middle row: same as the top, but now magnified 



y-axis, showing the MG triplet eigenvalue (in one color) that crosses the rising singlet eigenvalue 



of the MG subspace (in a different color). Bottom row: The MG eigenvalues rescaled so that their 
sum is unity. Largest 50 eigenvalues are shown. 



to these 5 eigenvalues along with identical vectors from the remaining N/2 particles form 
the N— particle MG subspace. 



Fig. [T] (middle row) shows the intersection of the largest eigenvalue in the MG (which is 
triply degenerate shown using one color) with the rising eigenvalues of the dimer MG sector 
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(shown using a different color). This is a robust feature for all N and an even number of 



spins in the subspaces. More of the eigenvalues corresponding to the MG subspace are seen 
in the bottom row which shows the rescaled eigenvalues in this sector. In the rescaled figure 



which shows only the MG eigenvalues, the most prominent ones are again few and the two 
that are shown correspond to a pair of triplets that seem to be coupled strongly. 

Fig. |2] shows the entanglement spectrum defined as — In (Ay) plotted against J 2 . This 



figure now highlights the small eigenvalues in the MG subspace and a clear separation is 
seen as those belonging to the MG subspace now decrease as the MG point is approached. 
There are several sharp peaks that are seen in these figures and their density increases as 
N does. 

These signal eigenvalues of the RDM that either go exactly to zero or come very close to 
it (it is sometimes difficult to tell with given numerical resolution in J 2) and interestingly 
resume their career immediately thereafter. 

It is then quite apparent that there are only few dominant eigenvalues of the density 
matrix, even away from the MG point. That these are actually those that produce the 
dimer is made clear by studying the entanglement spectrum of a superposition of the dimers 
\R) and \L). Towards this end consider the state 



\^mg) = oii\Rn) + a 2 \L N ) 



(5) 



where \Rn) and \L N ) are as defined in Eq. (T2J), and ai )2 are real. It is the simplest type of 
a "Valence Bond State", which is a superposition of dimerized states 421 ] . While in general 
such VBS states have been quite extensively studied, including from the point of view of 



entanglement 



21 



43M45|| to our knowledge a detailed analysis of the simple state in 



\^mg) at finite ./V and arbitrary partition sizes has not been reported. 

We begin here by evaluating the required RDMs. Let Na = 2k be the number of particles 
in the subsystem A (k is any appropriate integer > 1) whose density matrix is given by 
(details are relegated to an Appendix A): 



MG 

Pa 



ti B = al\R 2k ){R 2k \ + a\ 



Air \, j I _ hk 



+ 



1) 



(N-2k)/2 



\R2k)(L2k \ + \L2k)(R 



2k\ 



(6) 



where |i? 2 fc) 



2(A r -2fc)/2 

(1 2) • • • (2Jfe - 1 2k) and \L2k-2) = (2 3) • • • (2k - 2 2k - 1) while \L 2k ) 



\L2k-2) (2fc 1) are dimers of part A; \L2k-2) does not contain the singlet between the first 
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FIG. 2. The entanglement-spectrum that shows the MG as excited states (50 of the largest density 
matrix eigenvalues are plotted). The separation of the MG and MG eigenvalues is seen clearly 
here, as well as the proliferating number of eigenvalues that vanish at isolated points along J2 are 
seen as sharp peaks. 



and the "last" (2k) spin of part A. As the inner product {R K \L K ) = (-1) K/2 /2 K/2 ~ 1 , it is 
readily verified that tr(p^ /G ) = a\ + a| + 2a\a2(RN\LN) = mg\^ mg) ■ Thus if a\ and a<i 
are taken such that \^mg) is normalized the trace of the RDM pf G is indeed 1. 

To find the spectrum of p^ G , it is useful to express the identity operator in the space of 
spins 1 and 2k, Ii <g> I 2 k in terms of the complete set of corresponding Bell state projectors. 
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This results in: 



pT =al\R 2k )(R 2k \ + ^ \L 2k ){L 2k \ + ^ \£\L 2k )(L 2k \ j + 
_^_ ( _l)(^)/2 (li^X^I + \L 2k )(R 2k \) . 

Here \L l 2k ) are |L 2fc _ 2 > <g> |0<) 2M , where |#) = (|01) + |10))/V2, |0 6 2 ) = (|00) + |11»/V2, 
and = (|00) — \11))/V2 are three Bell states, the remaining one being the singlet that 
along with \L 2k _ 2 ) results in \L 2k ). These Bell-state augmented dimers are quite easily seen 
to satisfy the following: (R 2k \L l 2k ) = (L 2k \L l 2k ) = 0, for I = 1,2,3. Thus \L l 2k ) are three 
degenerate eigentates of p^ G with eigenvalues \ 2 MG = Af G = \f G = a 2 /4. Two other 
eigenstates are linear combinations of the nonorthogonal states \R 2k ) and \L 2k ), and the 
resultant eigenvalues are 



\MG _ ' 
Ai,5 "2 



(8) 

These 5 eigenvalues of the 2k— particle (k > 1) RDM are the only nonzero ones, and 
it is easily verified that they add up to the trace of the RDM. If the initial dimer state 
is normalized they add to unity. They are ordered according to their typical magnitude, 
especially when a\ = a\ = 1/2, with \^ G being the largest and \f G being the smallest 
eigenvalue. When N = oo, and a\ = a\ = 1/2, the eigenvalues are, as indicated earlier, 
(1/2, 1/8, 1/8, 1/8, 1/8) and the entropy or entanglement is 2ebits. For finite N the entropy 
is smaller, for example when N = 8 and 2k = 4, taking the normalized state \^mg) with 
«i = a 2 = 2/3 leads to \ 2 IG = \^ IG = \f G = 1/9 as the eigenvalues for the degenerate 
triple of states and the other two are Af G = (2 + >/3)/6 ~ 0.622 and Af G = (2 - V3)/6 w 
0.044, while the entanglement is 1.683 ebits. To take an example of an odd case, let 
iV = 10 and 2k = 4, the normalized state that has momentum n, maybe taken as \iPmg) — 
y/8/17{\R 10 ) - \L W )). Thus with a x = -a 2 = y/8/17, the above gives A^ 3 G 4 = 2/17 w 0.117 
and Af G = (11 + 2>/l9)/34 w 0.579 and Af G = (11 - 2 v / l9)/34 Pd 0.067. This may be 
compared with the eigenvalues of the RDM near the MG point. For example for iV = 10 
and 2k = 4 when J 2 = 0.4975, the triplets in the dimerized part have eigenvalues of 0.1175, 
while the large eigenvalue and the smallest one in this part are 0.582 and 0.065, which indeed 
compare well with the numbers derived above. 
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Thus the "entanglement-spectrum" at the MG point consists of only 5 levels. The fact 
;hat the ground state of the MG model for < J 2 < 1/2 does not undergo any crossings 
15| indicates a certain robustness that will be reflected in the entanglement-spectrum as 
well. As evidenced also by results shown in Figs. 1 and 2, indeed a 5-dimensional subspace 
dominates the entanglement spectrum and evolves to the one derived above when J 2 — >■ 1/2. 
Thus Eq. (j3J) maybe split into two parts with |^(J 2 )) = Y^=i V^K^jl^ji^)) A^ji^)) b + 
Ylj=t V ^ j ( ^2 ) I <\> j ( J 2 ) ) a 1 4> j ( J 2 ) ) b, which is the separation that is alluded to in Eq. (J3J). Thus 
\1pMGJ2)) is in the 5-dimensional "dimer" MG subspace that dominates the state, while 
IV'mg^)) belongs to the (2^ — 5)-dimensional subspace which constitutes the rest. The 
eigenvalues Aj(J 2 ) for 1 < i < 5 are defined as those that evolve to Xf 10 at the MG point 
J 2 = 1/2. Thus it follows that 

5 2 N A 

a\j 2 ) = ]T Xj(J 2 ), /3 2 ( h) = E \ S W = 1 - o?{ J 2 ), 



j=l j=6 



2 N A 



WmgW) = E \I*^\mj2))a\mj*))b, i^wgW) = E J^yi^^jxi^cJa)^. 



3=1 



(9) 



The identification of the eigenvalues belonging to MG is complicated slightly by the fact 
that the largest eigenvalue in this set crosses the eigenvalue that becomes X^ IG of Eq. (jS]). 
Indeed this "rising" eigenvalue in the dimer subspace is coupled to the largest eigenvalue 
state and its dominance in the spectrum seems correlated with the dimerization process. 

It is important to note that the identification of the iV-particle pure states \iPmg{J2)) 
and |V'mg(^2)) from the Schmidt vectors is dependent on the partition sizes Na and Nb and 



thus usage of terms like MG and MG subspaces, is predicated upon a definite partition 
dependence, usually the symmetric one, corresponding to = N B . 

Following the above considerations, one may find three entanglements between con- 
tiguous spins and the rest: 



2 N A 5 



s(j 2 ) = - E W2) MM J2)), SmoW = - E log (Wr 

i=i i=i v 1 x v 

„ m ¥^(J 2 ) fW 2 )\ 



(10) 



i=6 

whose interpretations respectively as the entanglements in the ground state (\^(J 2 ))), and 



separately in the MG and MG parts (IV'mg^)) and |V'mg(^2)}) of the ground state is 
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straightforward. The behaviors of (3 2 (J2) and S(J2) are shown in the Figs. ([3]) top left 
and right plots respectively, while the bottom left plot of Fig. (j3J) shows Sa/g^), the 
entanglement in the MG subspace projection. The bottom right plot shows Sjjq{J2). In 
these figures for even cases, (N/2 is even), Na = N/2, while in the odd case Na = N/2—1. It 
is interesting that while the entropies, S( J 2) and Smg are monotonic, the entropy Sjjq shows 
a minimum in the vicinity of the dimerization transition. If indeed these are entanglement 
signatures of this quantum phase transition, it is interesting that it is found in the "non-MG" 
part of the state. Of course this part increases in dominance as N increases, see Fig. (J3]) 
(top left plot). 

That there is a fairly significant dimerized part that is already present in the small N 
Heisenberg model maybe the reason why the reason why the entanglement signatures of the 
transition are not easy to see; but once the dimerized part is excised, at least in part, the 
remaining "grass" seems to reveal the transition. It should also be noted that calculations 
not presented show that if the entropy S is itself split into a MG and MG part without 
rescaling the eigenvalues, then these are monotonic on [0, 1/2]; the interpretation of Sjjq as 
a entanglement is necessary. 

It is also observed that the spectrum of the RDM for various partition sizes Na = 2k 
are qualitatively similar including the crossing of the lowest eigenvalue corresponding to the 
MG subspace with the triplet from the MG one. The existence of a minimum entropy of 
entanglement for the IV'mg^)) state between Na and the rest of the spins is interestingly 
a robust feature, as shown in Fig. (j3J). Of course for k = 1 (two-spin RDM) is special, the 
number of eigenvalues of the RDM being 4 and all remain non-zero at the MG point. In 
fact it is easy to see from the Eq. (]6]), which becomes a Werner state, that these eigenvalues 
correspond to the one dominant one and the triplet. The "rising" state is absent from this 
spectrum and is a property of chains with more than 4 spins. It is quite essential that the 
number of spins in the subsystems A and B are even. If there are an odd number of spins (in 
the subsystems) the number of eigenvalues in the RDM that are non-zero at the MG point is 
4 and the entropy Sjjq remains monotonic in [0, 1/2]. There is also the added complication 
that the ground state of the J\ — J 2 model in this range has zero momentum when N /2 is 
even and momentum ir otherwise. 

The eigenvalues in MG themselves have structure and a hierarchy that is not unlike that of 
the dimerized state. While the largest triplet MG eigenvalue, A 6 ( J 2 ) decreases monotonically 
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FIG. 3. The sum of the eigenvalues corrsponding to the MG subspace, /3 2 (J2), is shown on the 
top left plot for various values of the number of spins iV from N = 8 to N = 24 in steps of two, 
while the top right plot shows the entanglement of Na spins with the rest (the entropy S 1 ^)) in 
the complete ground state . The bottom left plot shows the entanglement of the projection in the 
MG subspace (the entropy Smg(^2))- The entropy Sjj^^), the entanglement of the projection 
in the MG subspace as defined in Eq. ()10p is shown in the bottom right plot. Na = N/2 for the 
even case and N/2 — 1 for the odd. 

in [0, 1/2], the scaled value (divided by /3 2 (J 2 ), see the bottom panel of Fig. ((T])) shows a 
single peak again in the vicinity of the dimerization transition, which maybe the origin of 
the minimum in the entropy of the grass. Indeed — log(Ag(J2)//3 2 (J2)) is the so-called min- 
entropy, S^jjq, and along with the von-Neumann entropy is a special case of the Renyi 
entropies. It is guaranteed from general considerations that S^j^ < Sjjq. 

To explore entanglement sharing in the pure iV— particle states |V>mg) aric ^ H'mg), for a 
fixed partition iV = Na + Nb, one can study multipartite measures and 2-spin measures such 
as concurrence. The latter is studied in the next section, while for the first, the entanglement 
entropy of m contiguous spins is numerically calculated, in the case of symmetric partitions 
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FIG. 4. The entropy Sj^q{J2) for N = 16 spins with different bipartition sizes Na — Nb that are 
indicated. This is then the entanglement of Na spins with the rest for the MG projection. The 
existence of a minimum is robust to altering paritition sizes. 

Na = Nb, the results of which are shown in Fig. (jSJ) for a chain of length N — 16. The 
translation symmetry is lost on projection and by construction, from the Schmidt decompo- 
sition, these states have a symmetry of shifting by N/2 sites. Thus it makes a difference as to 
where the first of the m contiguous spins is chosen. The case of "0 shift" corresponds to the 
first being also the first in the block of N/2 spins that remains after tracing. Further shifts 
refer to right shifting the first spin in the block by the indicated amount. The rather more 
complex entanglement sharing of the state \iPjyg) is seen here. The dependence on the shift 
of the first spin of the m blocks is clear and for no shift the prominent feature remains the 
m = 8 case that has already been discussed above. However the case of m = 4 also shows a 
local minimum, albeit a shallow one, in the entanglement at exactly the same value of J2 as 
for m — 8. The other values of m also indicate the fair amount of multipartite entanglement 
present in this state. Shifting the first spin away, now explores different entanglement fea- 
tures, for instance with a shift of 1, and m — 2, this is the entanglement of spins 2 and 3 with 
the rest when the whole state is IV'mg)- The first observation is that there is no minimum 
anymore for any value of m, especially 8 spins. Thus, it is required for the minimum in 
the entropy that the block coincides with the partitioning in the Schmidt decomposition. 
The second is the considerably large entanglements that is present, for example with blocks 
shifted by 3. These are in sharp contrast to the behavior of the corresponding quantities 
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0.2 0.3 
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FIG. 5. The entanglement entropy between m contiguous spins and the rest for the state \tPjjq) 
Here N = 16, and the "shift" refers to the first spin site in the block of m spins. 



for the state \iPmg)i wherein the entanglements are seen to be monotonically increasing and 
the shifts do not change the features much, and the entanglement entropies are only about 
half as large. 



III. CONCURRENCES IN THE SCHMIDT VECTORS AND THE MG, MG 

PROJECTIONS. 

Attention is now turned to a more detailed study of two-spin entanglements. In particular 
one wishes to know the nature of entanglement in the eigenstates of the RDM of p^ A . 
These correspond to the A^-particle Schmidt vectors in a Schmidt decomposition of the 
ground state. Also of interest is the concurrences present in the corresponding projected 
N— particle pure states \iPmg) and The concurrence is a one-to-one function of 

the entanglement of formation and the recipe to obtain the concurrence between any pair 
of spins which are either in a pure or a mixed state p is as follows [24J: compute the 
eigenvalues of the matrix p[a y (g) a y )p*(a y <g> cr y ) (the complex conjugation being done in 
the computational basis). The eigenvalues are guaranteed to be positive and if they are 
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arranged as {Ai > A 2 > A 3 > A4}, the concurrence between the pair of spins considered is 
given by C = max[0, y/Xi — \f\2 — y/Xs — x/Al] . The concurrence C is such that < C < 1, 
with zero for the case of an unentangled state and unity when it is maximally entangled. 

Recall that the MG subspace for a given even partition is spanned by 5 states whose cor- 
responding eigenvalues are a large and decreasing one, three degenerate and nearly constant 
ones, while the last is small and increasing. Concentrating on the case Na = N/2 and N an 
even integer, it is sufficient to study the N/2— particle pure states \4>i(J2))a, 102,3,4 (^2)) a, 
and \4>§{J2))a ( see Eq. ([9])) respectively. Collectively they contribute to the normalized state 
\iPmg{J2))- The complementary subspace is the normalized state l^gc^)) whose principal 
contribution comes from a triplet whose eigenvalue is decreasing and intersects with the 
increasing lowest eigenvalue from the MG subspace. In all of these states one can look at 
the nature of pairwise entanglement via nearest neighbor pairwise concurrence j^j], which 
is a genuine and well-used measure of entanglement between two qubits or spin-1/2 parti- 
cles, especially useful when they are not in a pure state. One may study the concurrences 
C(i,i+i) (J2) between spin at i and i + 1 (identifying L + 1 as the first spin), as well as their 
totals either over the entire chain, or over two parts, where i is even or when it is odd. 
Note that while the ground state has translational invariance, this is typically broken in the 
states \i[>mg(J2)) and |V ; mg(^2)), and the corresponding Schmidt vectors. Thus (7(^+1) ( J 2) 
are typically different for different values of i, unlike in the original state. 

As the dimer part of the state survives till the MG point, it is likely to have large pairwise 
concurrences. The top panel of Fig. ([6]) shows these for the most dominant state, namely 
\4>i(J2))a and the rising state |05(J 2 ))a- It is seen that the nearest neighbor concurrences 
show a clear progress to dimerization as J2 increases. In the case of |0i(J 2 )),4, the entan- 
glement between the alternate bonds starting from the first is large and increases with J2, 
while the others decrease and vanish well before the MG point, and in this respect is like 
\Rn/2)', while for \4>b(J2))a the highly entangled bonds start from the second spin, and in 
this respect is like \L^/ 2 ). The insets show the sum of the concurrences in the even and 
odd sub lattices of the N/2 spin chain, and it is seen that the entanglements vanish in the 
alternate bonds again in the vicinity of the dimerization transition, but not at exactly one 
point. Also notice that for |</>s(J2))a entanglement develops for distant spins at sites 1 and 
8 for the case shown of N = 16. The full chain has been "cut" keeping sites 1-8 and trac- 
ing out 9-16. This singles out the sites 1 and 8; also note the reflection symmetry that is 
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N=16 - for the state I (pfi 2 )) A N=16 - for the state I <P 5 (J 2 )> A 

~3 * I , sum ol pair concurrences' 




j 2 j 2 

FIG. 6. The concurrence between neighboring spins in the Schmidt vectors, which are pure states 
of N/2 spins, for the symmetric partition case of N = 16. The top left is for the vector |0i(Ja)) 
corresponding to the highest eigenvalue, top right is for the vector |0s( corresponding to the 
rising eigenvalue, while the bottom left is for the mixed state P234 corresponding to the triplet. 
The bottom right plot is for the mixed state P678 corresponding to the largest triplet of eigenvalues 



corresponding to states in the MG subspace. The insets in the corresponding graphs show the sum 
of the alternate pair concurrences for each of the above described eigenstates. 



apparent from the distribution of concurrences amongst the spins 1-8. It maybe noted that 
the concurrences for the case of the ground state do not show such structures that reveal 



the dimerization 



15 



26 



27| 



In case of the triplets (as in the MG subspace, or the states corresponding to the largest 



eigenvalue in the MG subspace) the states are not unique due to degeneracy. However 
the projector onto the degenerate three-dimensional subspace is unique and can be used to 
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define a density matrix, for example for the triplet in the MG subspace consider the state: 

P234(J 2 ) = \{\HJ2))A A (<hW\ + \HJ2))A AifoWl + \MJ2)) A A (M J *) I ) > ( U ) 

and the corresponding mixed state for the most prominent triple of MG, say P678- The 
concurrence in the bonds of these states are shown in the lower panel of Fig. (J6]). The state 
p 234 displays large entanglements in the alternative bonds starting from the second spin and 
is in this respect like the rising state | <^>5 ( ^2 ) ) , except that there is here no entanglement 
between the distant spins 1 and 8. This distinction from the rising state is understood on 
calculating the two-qubit reduced density matrix of the triplet density matrix at J 2 = 1/2: 

(p 23 4(J 2 = 1/2)) 18 = i (|00)(00| + |11)(11| + I + +)(+ + I) (12) 

the reduced density matrix of the spins 1 and 8 is (p 2 34(</2 = l/2))is, the state |+) = 
(|0) + |l))/y / 2- The separability of these spins in this density matrix is then evident, and 
entanglement appears to be absent not only at J 2 = 1/2 but for the entire range considered. 
The entanglement of other bonds decrease from the Heisenberg point and vanish again well 
before the MG point. Indeed the point where the dimer states takes on a pure alternate 
bond entanglement is again in the region of the dimerization transition. 

A similar analysis for p678 is shown in the same figure and presents a somewhat different 
picture, with the dimerization not being uniformly present. While the concurrence between 
3 — 4 and the symmetric 5 — 6 spins are large, the rest of the nearest neighbor entanglements 
are nearly zero. The concurrences that is present in the 3 — 4 pair is also decreasing from 
the Heisenberg chain as the MG point is approached. It is observed that the entanglement 
between the 1 — 2 and 7 — 8 pairs, which starts at zero, develops as J 2 increases and is non- 
zero at J 2 = 1/2. This indicates the existence of some dimerization in the MG subspace, but 
of a different kind than in the MG subspace. The effective overall decrease of the two-qubit 
entanglements is in sharp contrast to that found for the Schmidt vectors that span the MG 
subspace. 

To analyze this further, pair concurrences were calculated for the N spin state [ipMai^)) 
and the nearest neighbor pair concurrence is shown in Fig. (J7|). While this looks similar to 
the case of the state P678, which is indeed the dominant part of the prominent 

difference is the somewhat large entanglement between the 1 — 16 and 8 — 9 pairs (here 
N = 16) which exists for the Heisenberg chain and which decreases away as the MG point is 
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FIG. 7. The graph shows the nearest neighbor concurrence in the N = 16 spin state ipjjci-h)- 
The spin pairs are indicated above their corresponding curves and the inset is a magnification. The 
pairs not indicated, namely (1, 2), (7, 8), (9, 10), (15, 16), have zero concurrence throughout. 



approached. Once again the overall decrease in the concurrence is in contrast to that for the 
MG subspace and is consistent with an increase of the von Neumann entropy in as much 
as one can think of monogamy of entanglement being operative and the entanglement be- 
comes of a more multipartite kind. Indeed the structure of even the two-spin entanglements 
present in this state is neither of a \L^) kind nor of a \Rn) kind, but rather a mixture with 
some bonds being either very weakly or not at all entangled. For example the entangle- 
ment present in 1 — 16 and 8 — 9 is consistent with a \Ln) kind of dimerization, while the 
prominent entanglement between 3 — 4, and 5 — 6 resembles \Rn)- It is interesting that as 
the dimerization progresses, even in this "grass" contribution there is a tendency to choose 
a type of dimerization, with the |Ljy) kind taking a backseat at around the dimerization 
transition. However also the decrease in the concurrence is consistent with a rising entropy 
as observed in Fig. (|3]) (bottom right plot). The principal features discussed above have 
been verified to remain intact for the case of N = 24 spins. 
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IV. DISCUSSIONS AND SUMMARY 

In this paper the frustrated J\ — J 2 antiferromagnetic Majumdar- Ghosh model has been 
revisited with a view on entanglement properties, both multipartite and those between pairs 
of spins. Entanglement studies of the ground state that have revealed signatures of the 
dimerization transition have hitherto relied on scaling of the entropy with the system size. 
However in this paper several suggestive simple signatures are presented, from those that 
involve von-Neumann entropy to concurrence between spins. For this the principal tool is 
the well-known Schmidt decomposition that combined with the existence of the unique MG 
point (J2/J1 = 1/2) provided an opportunity for a projection of the ground state into two 
orthogonal subspaces that are unique once the partition in the Schmidt decomposition is 
fixed. The dominant subspace is only 5 dimensional and contains the complete state at 
the MG point. The complementary subspace whose significance wanes from the Heisenberg 
point (J2/J1 = 0) contributes only marginally, but this contribution increases with N, the 
number of spins; for instance for N = 16 this contribution is roughly 3% when J2/J1 = 0. 
These subspaces are indicated as MG and MG, although again the partition dependence in 
the Schmidt decomposition is implicit. 

The entanglement between spins and the rest of the spin chain is known to have 
different scaling laws as criticality is lost at J 2 j J\ ~ 0.24. What is shown above is that while 
the entanglement of the full or dominant projection in the MG subspace is monotonic, the 
projection onto MG has a minimum in the vicinity of the transition for even Na, at least 
for the values of N that has been explored. This feature is robust against various different 
partitions of the ground state. How robust this feature is to increasing number of spins 
remains to be seen. 

To understand better the behavior of entropies, one may calculate further the entan- 
glements present after the states \L N ) or \Rn) are projected out from the ground state 
|*(J 2 ))- That is the quantities S L (J 2 ) = 5((|*(J 2 )) - \R N )(R N \^(J 2 )))t) and S R (J 2 ) = 
S((\^f( J 2 )) — \L N ) (L N \^/(J 2 )))t), the entanglements of N/2 contiguous spins with the rest 
(t is a normalization constant), are found and plotted in Fig. ((H]). Note that Sl(J 2 ) tends to 
the entanglement of the state \Ln) at the MG point, namely 2, while Sr(J 2 ) tends to that 
of \Rn), namely 0. It is interesting to note that the entanglement of these symmetry broken 
states are now monotonically decreasing already for small values of N, reflecting well the 
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FIG. 8. The etanglements of half the chain with the rest for Sr{J2) and Sx(J2) after projecting 
out the fully dimerized \Ljq) or \Rn) states respectively. Shown also is the case for the complete 
ground state entanglement S^J^) for comparison, and here N = 16 spins are used. 

fact that the entanglement sharing in the spins is changing from a more complex situation at 
the Heisenberg point to the fully dimerized one at the MG point. The dimerization, which 
is leading to the formation of couples that are unentangled with any other spin discourages 
multipartite entanglement. 

An important complementary view of entanglement sharing is provided by calculating 
the concurrence between pairs of spins. The state in which these are measured are however 
not the ground state itself, but the eigenf unctions of the reduced density matrix, or the 
Schmidt vectors. The Schmidt vectors (now pure states of Na spins) of the MG subspace 
show a clear progress towards dimerization as J<i increases with the the most dominant 
state resembling \Rn/2) an d the rising state, \L^/ 2 )- The triplet also shows dimerization 
as in |Ljv/2) except for the end spins being unentangled. Here "dimerization" is seen as 
the vanishing of concurrence on a sublattice, while the complementary one develops into 
pairs with maximum concurrence. The most dominant eigenvalue corresponding to triply 
degenerate states in the MG subspace was also studied using pair concurrences and it 
presents a different picture compared to the 5 Schmidt vectors in the MG subspace, in that 
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the concurrences tend to decrease as the MG point is approached. The projection of the 
state on the iV particle MG subspace, IV'mg^)) also shows interesting differences and larger 
multipartite entanglements. The initial decrease of the entropy Sj^q which contributes to 
the non-monotonic character of this entropy may have its origins in the overall tendency for 
decreasing entropy as evidenced on projecting out the |Ljv) or \Rn) states, however more 
study is warranted on the exact origins and significance, if any, of this. 

If a bipartite split with one block containing the spins at odd sites and the other block 
containing the spins at even sites ("comb entanglement") is taken, it presents a complex 
entanglement spectrum with many crossings, and while this is interesting, the dimerization 
transition seems difficult to unravel. Also the present study has calculated non-nearest 
neighbor concurrence in the various states presented, but most of them are indeed zero. 
Preliminary investigations of the J\ — J2 model with quenched disorder in J 2 reveals a 
certain robustness of the above analysis. Small disorders lead to the exact crossing at the 
MG point being replaced by an avoided one, and there is still to a large extent only a 5 
dimensional dominant space and hence a split into an MG and MG subspaces persists. 
Further study is needed on how non-dimerized subspaces such as MG dominate in the large 
N limit. 
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Appendix A: Computation of the reduced density matrix and its spectrum at the 

MG point 

There is more than one way to derive the reduced density matrix, and in the following 
a direct approach is used. First start with the superposition \ip) = oi\Rn) + ch2\Ln) of the 
two dimer states that are eigenfunctions at the MG point, namely \Rn) and |Ljv) (as given 
in Eq. ©). 

On taking a bipartite split of 2k contiguous spins, the block A, and the remaining N — 2k 
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spins, the block B. One traces over the block B to find the reduced density matrix: 

P2k = tT B (\ip)(ip\) = ajp 2kl +alp 2k2 +a 1 a 2 (p2k- i + pL 3 ), ( A1 ) 

where p 2kl = tr fl (|J?^)(i?jv|), p 2fea = tx B (\L N )(L N \), and p 2fc3 = tr jB (|i? A r)(L Af |). 

Denote now \R 2k ) = (12)(34) • • • (2k-l 2k) and |L 2fc „ 2 ) = (23)(45) • • • (2k-2 2k-l); these 
being spin states which are not affected by the partial tracing operation. It is straightforward 
to calculate p 2kl as no singlet "bonds" are cut due to the structure of \Rn)- However for 
calculating p 2k2 the singlets between (Nl) and (2k 2k + 1) are broken which results in 
maximally mixed states — and — at the ends. The remaining tensor products of singlets 
\L2k-2) are left unaffected. Thus it follows that 



a 



P2fcj = a l \R 2 k)(R2k\, P2k 2 = "^y 1 ® \ L 2k-2}{L 2k -2\ <8> hkj ■ (A2) 

The remaining part involves cross terms, which is written explicitly by introducing standard 
o z basis for the spins in block B: 

P2k 3 = (i2k+\- ■ -iN\RN)(L N \i 2k+ i- ■ -i N ) (A3) 

isk+i,— ,*jvS{0,1} 

It is easy to verify that (iii 2 |(|01) — |10)) = (— l) 1 5i 1 ,i a ©i i where the © denotes an addition 
modulo 2. 

(— 1) *2fc+l+»2fc+3H NjV-1 

and a similar expression is found for (LN\i 2k+ \i 2k+2 ■ ■ - In) which involves the untraced part 
\L 2k - 2 ) as follows, 

/ 1 v ((AT-2fc)/2)-2 

(-kjVl*2fc+l*2&+2 ■ ■ • iJV/ = (-^2*— 2 ] ] <^2fc+2i2fc+3©l ' ' ' ^tjyr_ 3 ijv-i0l X 

((2A;2A; + l)|i 2fc+ a)<(JVl)|i JV > (A5) 

The "end spins" are taken into account as ((2k 2k + l)\i 2k+1 ) = ^((0| 2 fc<5i ii2fc+1 - (l\ 2k 5 0:i2k+1 ) 
and ((iVl)|zjv) = -j%{(l\i8o,i N — (0 jx^i,ajv) - Using this expression along with Eqs. (1A4|A5j) 
and substituting them in Eq. (IA3j) . the final form of p 2k (after some straightforward algebra 
taking care of the modulo 2 addition) is found to be 



p 2k = a\\R 2k )(R 2k \ + a\ 



^1 „ I r \l T I „ ^2k 

— <8> \L 2k - 2 )(L 2k ^ 2 \ g) — 



Ql^2 / ,x(JV-2fc)/2 



|-R2fc)(-^2fc| + |-^2fe) (-R2fc| 



(A6) 
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It is to be noted from the structure of p 2 & that the coherent term is only of the order of 
2~\ N ~ 2k >/ 2 and hence exponentially decreases with the number of spins in block B. 

The eigenvalues of p2k are now calculated. For notational convenience, for K spins, 
define p = (R K \L K ) = (_l)*/2 2 (i-K/ 2 ) an d 7 = 1 2 V_ 2fc)/2 • Writing the h ® hk 



in lli <8> \L2k-2) (L2k-2\ <8> hkj as the sum of projectors into the four Bell basis, \<j> l b ) (1 < 
I < 4) we can define new states \L l 2k ) = |L 2 fc- 2 ) ® \<t>b)^k,u nere ^ = 1)2,3,4, we then use 
the (easily obtained) properties that {R2 k \L l 2k ) = (L2k\L l 2k ) = 0, I = 1,2,3. Explicitly: 
\L\ k ) = |L 2fe _ 2 )((|01) + |10})/V2) 2fe)1 , |LL> = |L 2fc -2>((|00> + |H))/V2) 2M , \L\ k ) = 
\L2k-2) ((|00) - |11)W2) and = |L 2fc _ 2 ) ((|01) - 10))/V^) = |L 2fc ). Now using 



the above we can rewrite p2k as 

3 

'+7 



P2k — Ol{\R2k) (R2k\ + —r\L2k) (L 2 fc | + — p 



4 



|-^2fc)(L 2 A;| 



1=1 



\R2k){L 2 k \ + |L 2 fc)(_R 2 fc| 



(A7) 



It is straightforward to verify that |L 2fe ),|L 2fe ) and |L 2fc ) are three eigenstates of p2k with 
degenerate eigenvalues This corresponds to the triply degenerate eigenvalue in the 

entanglement spectrum at the MG point. It is clear that the other eigenvalues correspond 
to eigenvectors in the two dimensional subspace spanned by \R2kl and \L 2 k), and hence there 
are only two of these. Either defining the orthogonal vectors |-R 2 fc) ± \L2kl 01 proceeding 
to define linear superpositions of \R-2k) and \L 2 k) as the eigenvectors, a straightforward (if 
somewhat lengthy) calculation leads to Eq. ([H]). These two eigenvalues correspond to the 
most dominant state and the rising state in the entanglement spectrum of p 2 ^. 
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